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Abstract 


A  brief  Introduction  to  continuous  mathematical  control  theory  is 
presented  first.  Then  a  model  havi’'-  two  state  variable  accounts,  cash  and 
wheat,  is  defined  by  means  of  differential  equations.  Adjoint  functions, 
the  Hamiltonian  function,  and  the  optimum  (bang-bang)  policy  are  derived. 
Interpretations  of  the  Hamiltonian  and  switching  functions  and  a  numerical 
example  are  then  given.  Finally  the  results  are  extended  to  the  case  where 
transaction  costs  and  spoilage  costs  are  incurred. 
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APPLICATIONS  OF  MATHEMATICAL  CONTROL  THEORY  TO 
ACCOUNTING  AND  BUDGETING 
I.  THE  CONTINUOUS  WHEAT  TRADING  MODEL 


Yuji  Ijiri  and  Gerald  L.  Thompson 


I.  INTRODUCTION 

In  this  paper  we  shall  apply  mathematical  control  theory  to  the 
analysis  of  a  simple  model  of  an  accounting  and  budgeting  problem.  Our  ob¬ 
jective  is  not  only  to  make  the  model  operational  for  practical  applications 
but  also  to  use  the  framework  and  results  of  mathematical  control  theory  in 
order  to  get  new  insights  into  accounting  and  budgeting  problems.  In  the 
present  paper  we  shall,  however,  limit  our  objective  to  the  latter  point, 
leaving  the  former  objective  to  our  follow-up  paper  [3]  which  deals  with  the 
discrete  wheat  trading  model. 

We  shall  be  primarily  concerned  here  with  the  process  of  the  operations 
of  the  firm  looked  at  from  an  accounting  standpoint.  In  accounting,  a  firm 
is  represented  by  a  set  of  quantities  that  indicate  the  stocks  of  assets  of 
various  types  that  are  under  the  control  of  the  firm,  and  the  firm's 
activities  are  indicated  by  the  changes  in  these  quantities.  Using  mathe¬ 
matical  notation  we  may  say  that  the  firm  at  time  t  has  a  vector  x(t) 
whose  components  indicate  the  physical  amounts  of  each  asset  the  firm  has  at 
time  t.  The  results  of  the  activities  of  the  firm  between  tj  and  t^  •’re 
measured  by  the  vector  difference  x(t£)  -  x(t^).  From  a  conceptual  point 
of  view  it  is  convenient  to  consider  derivative  x(t)  defined  as, 

K<t)  -  -jr^  -  r  t;  :  t, —  • 
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Then  we  can  say  that  the  firm's  activities  dt  time  t  are  indicated  by  x(t), 
which  in  turn  affects  x(t). 

The  activities  x(t)  of  a  firm  are  aimed  at  changing  the  state  x(t) 
of  the  firm;  but  thev  are  constrained  by  the  state  of  the  firm.  For  instance, 
the  production  rate  at  time  t  may  be  constrained  by  plant  capacity  at  time 
t;  interest  paid  out  depends  upon  the  amount  of  loans  outstanding,  etc.  Such 
relationships  between  assets  x(t)  and  activities  x(t)  need  not  be  fixed. 
They  may  be  changed  exogenously  by  factors  outside  the  control  of  the  firm, 
or  they  may  be  changed  as  the  result  of  th.»  intentional  effort  of  the  firm. 

For  instance,  plant  deterioration  occu.3  1  ecause  of  weather,  machine  wear 
and  other  natural  causes,  but  it  may  be  corrected  by  the  firms  engaging  in 
repair  activity.  The  prices  of  raw  materials  may  be  determined  exogenously, 
but  the  firm's  sale  price  may  be  partially  or  wholly  at  Its  own  discretion. 

When  there  is  freedom  in  the  asset-ai  rvity  relationship  that  the  manage¬ 
ment  of  the  firm  may  exercise,  it  implies  that  management  has  responsibility 
to  exercise  the  freedom  in  the  most  favorable  way  for  the  firm.  It  must 
answer  questions  such  as:  Should  the  plant  b»  operated  in  two  shifts?  Should 
the  plant  be  repaired?  Or  should  the  wareh  :se  be  modernized? 

In  answering  these  questions,  it  is  necessary  to  take  into  account  not 
only  the  production  side  of  the  problem,  but  simultaneously  the  marketing 
side,  the  finance  side,  the  personnel  side,  etc.  The  accounting  system  offers 
a  framework  by  which  such  aspects  of  the  firm  and  its  activities  may  be 
simultaneously  considered.  This  point  has  already  been  demonstrated  in  the 
paper  by  JJiri,  Levy,  and  Lyon  [2],  in  which  linear  programming  techniques  are 
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used  In  budgeting  and  financial  planning  based  on  the  double-e.  try  bookkeeping 
system. 

The  present  paper  is  similar  to  the  above  paper  [2]  in  approach  and 
objective,  however,  the  mathematical  technique  we  shall  employ  here  is  that 
matheu«ii.icai  concrol  theory.  In  the  next  section  we  s^tll  "ive  a  very 
brief  introduction  to  this  theory,  but  only  enough  to  indicate  the  steps  we 
are  going  to  take.  A  reader  who  wishes  to  go  into  the  reasons  behind  these 
steps  should  consult  some  of  the  references  to  mathematical  control  theory 
that  appear  in  the  bibliography.  See  [1,  4,  5], 

2.  AN  OUTLINE  OF  THE  MAXIMUM  PRINCIPLE 

We  shall  give  a  brief  outline  of  the  continuous  maximum  principle,  in 
a  form  general  enough  for  applications  to  the  rest  of  the  paper.  Consider 
the  problem 

P  T 

(1)  Maximize  f  F(x,m,T)dT 

x,m  0 

subject  to  the  differential  equation  constrains 

(2)  Xj,  =  f  ^ (x,m,  t)  ,  x^O)  =  x° 

for  i  =  l,...,n.  Here  the  n-component  (column)  vector  x  consists  of  the 
state  variables ,  and  m  is  a  k-component  vector  of  control  variables.  It  is 
assumed  that  m  belongs  to  some  prespecified  set  M  of  admissible  controls, 
and  we  usually  assume  that  M  has  certain  properties,  e.g.,  convexity.  The 

Q 

functions  F,  f,  and  the  inifial  values  x  are  assumed  known.  The  problem 
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is  to  choose  a  sequence  of  admissible  control  variables  m  from  M  in  order 
that  the  resulting  trajectory  in  the  state  space,  determined  by  (2), 
u.aximizes  the  objective  function  (1). 

We  now  proceed  to  put  the  problem  in  a  certain  standard  form.  Define 
a  new  variable, 

(3)  xii+1  =  JC  F(x,m,T)dT,  x^t'O)  =  °* 

Then  it  is  clear  that 

(4)  =  F(x,m,t)  =  f^fx.m.t),  x  ( 0 )  =  0  = 

This  defines  the  function  f  , ,  and  the  value  x°,,.  We  now  can  restate 

n+1  n+1 

the  problem  in  (1)  and  (2)  as 

(5)  Maximize  6x  =  x^+j(T)  *  6  =  (0,..,,0,1) 

subject  to 

(6)  =  f^(x.m.t)  ,  x1(0)  =  x° 

for  i  =  l,...,n,n+l.  This  is  the  Mayer  form  of  the  problem.  Note  that  the 
objective  function  is  now  a  linear  function  of  the  state  variables. 

We  now  introduce  n+1  adjoint  variables  (functions)  and  require 

that  they  satisfy 

n+1 

(7)  zx  =  E  z  x  =  k, 

i=l 

where  k  is  a  constant,  for  all  values  of  t.  Differentiating  (/''  with 
respect  to  time  we  obtain 
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n+1 

(8)  zx+zx  =  T,  (z.x  +  z  x, )  =  0 

1=1  11  1  1 

We  now  define  the  Hamiltonian  function  as 


(9)  5>=f"(x,z,m,t)  =  -zx  =  zx  =  zf 

where  the  second  equality  follows  from  (8)  and  the  third  equality  follows 
from  (6).  Differentiating  (9)  with  respect  to  x^ ,  yields 


(10) 


The  boundary  conditions  for  the  adjoint  variables  will  appear  later.  Note 
also  that,  differentiating  (9)  with  respect  to  z^  we  obtain 
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Note  that  condition  (7)  is  a  kind  of  "perpendicularity"  condition 
connecting  the  state  variables  and  the  adjoint  variables.  In  many  ways  the 
adjoint  variables  play  the  same  kind  of  role  in  control  theory  that  the  dual 
variables  play  in  linear  programming.  However,  adjoint  variables,  being 
functions  over  time,  are  inherently  more  complicated  than  dual  variables. 

We  are  now  in  a  position  to  state  the  Pontryagin  Maximum  Principle, 
as  follows:  In  order  to  solve  the  problem 

(12)  Maximize  Sx 
subject  to 

(13)  x  =  ,  x(0)  - 

1  '•  ~  l 


a  necessary  condition  is  that  there  exist  adjoint  functions  z  satisfying, 


(14) 


Maximize  -zf  — 


The  reader  ia  not  expected  to  accept  this  result  as  obvious,  and  we  are  not 
going  to  prr-,’,c2  a  proof  of  the  result  here.  However,  we  are  going  to  use  the 
above  procedure  to  solve  a  management  science  problem  which  we  call  the  con¬ 
tinuous  wheat  trading  model. 

There  are  several  iiuportant  things  to  note  about  this  theorem.  First 
of  all  it  gives  necessary  conditions  only.  Second  it  does  not  provide  a 
computational  method  for  finding  the  adjoint  functions  --  it  merely  gives  a 
way  of  checking  when  they  might  be  the  correct  ones.  Third,  notice  that  (13) 
is  an  initial  value  problem  for  the  state  variables,  while  (15)  is  a  terminal 
value  problem  for  the  adjoint  variables.  Since  we  must  have  a  simultaneous 
solution  to  both  problems  we  have  what  is  called  in  the  literature  a  "two- 
point  boundary  value  problem."  Such  problems  have  been  found  difficult  to 
solve  in  general  Fortunately,  in  our  case  ve  shall  be  able  to  provide 
solutions. 


5.  THE  CONTINUOUS  WHEAT  TRADING  MODEL 

In  order  to  obtain  insight  into  hew  mathematical  control  theory  can  be 
used  in  accounting,  let  us  consider  the  following  example.  A  firm  engages 
in  buying  and  selling  wheat.  The  firm's  assets  are  of  only  these  two  types, 
cash  and  wheat,  whose  balances  at  time  t  are  represented  by  c ( t >  and 
w(t),  respectively,  where  c(tl  is  measured  in  dollars  and  v(t)  In  bushels 
of  wheat.  The  price  of  wheat  from  time  0  to  lime  T  is  assumed  to  be  known  in 
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advance  with  certainty  and  is  given  by  a  function  p(t).  We  shall  let 
P  =  p(T)  be  the  terminal  price  of  wheat.  We  assume  that  the  Initial  afisets 
c(Q).  =  Cq  and  wvO)  =  W^  are  kncwn  constants. 

The  firm's  objective  is  to  buy  and  sell  wheat  during  the  time  0  to  T 
so  as  to  maximize  the  total  value  of  assets  at  time  T,  the  end  of  the  planning 
period.  Hence  it  always  values  its  current  wheat  inventory  with  the  terminal 
price  P,  in  making  the  trading  decisions.  Thus  the  total  value  of  the  firm's 
assets  at  time  t,  denoted  by  x(t),  is  given  by 

( 16)  x(t)=c(t)+Pw(t), 

We  denote  the  rate  of  buying  or  selling  wheat  by  m(t),  expressed  in  bushels 
of  wheat  purchased  if  positive  (or  sold,  if  negative)  per  unit  of  time.  We 
assume  that  m(t)  is  constrained  by 

(17)  M  ^  m(t)  £  M  for  0  S  t  s.  T, 

where  M  <  G  and  M  >  0  are  given  constants.  As  long  as  this  constraint  is 
satisfied,  the  firm  is  permitted  to  be  either  long  or  short  in  either  cash  or 
wheat,  or  both.  A  long  position  is  indicated  by  a  positive  stock  of  an  asset 
and  a  short  position  by  a  negative  stock. 

There  are  three  factors  that  change  the  cash  balance.  One  is  the 
Interest  received  on  cash  (or  paid  on  loan)  balance  at  a  given  constant  rate 
;;  per  unit  period,  which  is  to  be  accrued  continuously  and  to  be  received 
tor  paid)  in  cash.  That  is,  the  interest  on  cash  or  loan  balance  is  com¬ 
pounded  continuously  at  rate  R.  The  second  factor  is  the  inventory  holding 
cost  h(v)  that  must  be  paid  on  an  inventory  of  size  w(0  at  time  t.  The 
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function  h  is  defined  for  both  positive  and  negative  w  and  is  assumed  to 
be  an  increasing  function  of  )w|  .  Thus,  we  interpret  the  inventory  holding 
cost  to  be  the  cost  of  being  short  in  wheat  when  w(t)  <  0.  The  third  factor 
is  the  cos t  (or  revenue)  p(t)m(t)  of  purchases  (or  sales)  of  wheat  at  time 
t.  Thus  the  differential  equation  that  expresses  the  manner  in  which  cash 
changes  is 

fir 

(18)  c(t)  =  Rc(t)  -  h(w)  -  p(t)m(t). 

Contrary  to  the  cash  balance,  there  is  only  one.  factor  that  changes 
the  stock  of  wheat,  namely  the  purchase  or  sale  of  wheat,  m(t)  at  time  t. 

Thus  the  differential  equation  tor  w  is 

(19)  w(t)  =  m(t) . 

To  simplify  the  notation,  ,.e  shall  drop  (t)  from  all  expressions  with 
the  understanding  that  all  lower  case  letters  (except  e)  represent  functions 
of  time  unless  other  independent  variables  are  specified  (e.g.,  h(w)),  and 
that  all  upper  case  letters  represei.  constants. 

To  put  this  problem  into  standard  form  we  differentiate  (16)  obtaining 

(20)  x  =  £  4-  Pw 

=  Rc  -  h(w)  -  pm  +  Pm  * 

The  latter  expression  was  obtained  by  using  (18)  and  (19) 

Then  the  firm's  objective  is  expressed  by 

(21)  Max  *  c(T)  4-  Pw(T) 

subject  to  differential  equations  (18),  (19),  and  ^20)  and  constraint  (17). 

Tbs  system  of  eoixditions  (l/)-(21)  is  called  the  primal  system.  The 
variables  c,  w,  and  x  are  state  variables ,  and  ro  is  the  control  variable. 


9 


We  want  to  choose  a  sequence  of  values  for  the  control  variable  over  time  so 
as  to  maximize  the  terminal  value  of  the  firm's  assets. 

Following  the  outline  of  the  previous  sect  ion  we  introduce  ad  joint 
variables  for  each  constraining  equation.  Specifically,  we  choose  variable 
^r(t)  for  (18),  0(t )  for  (19),  and  the  constant  function  -1  for  (20).  [If  we 
had  chosen  an  arbitrary  function  for  (20)  it  is  easy  to  show  that  it  must  be 
equal  to  -1  everywhere;  we  omit  the  proof  of  this.]  Then  the  Hamiltonian 
function  is 


(22)  =  ijr  (Rc  -  h(w)  -  pm)  +  0m  -  (Rc  -  h(w)  -  pm  +  Pml 
=  (ij;  -  1)  (Rc  -  h(w))  -  [P  +  p(^f-l)  -  0jra 

We  shall  call  the  coefficient  of  m,  namely, 

(23)  s  =  P  +  p(*-l)  -  0 
the  switching  function. 

The  maximum  principle,  previously  discussed,  provides  as  a  necessary  con¬ 
dition  for  an  optimal  policy  that  the  control  variable  m  should  be  chosen  so 
es  to  maximize  Since  m  occurs  linearly  in  it  follows  that  the 

optimal  policy  is  bang-bang  and  is  of  the  following  form: 

f- 


(24) 


-  4 


M  if  s  >  0 
undetermined  if  s  =  0 
M  if  s  <  0 


The  case  where  s  —  0  is  called  s ingular  control.  For  most  problems  this 
occurs  infrequently  and  can  usually  be  ignored.  We  will  give  examples  later. 

Although  the  maximum  principle  gives  only  necessary  conditions,  the 
policy  given  in  (24)  is  unique.  Hence  there  is  only  one  extremal  policy  which 
therefore  must  be  the  optimal  policy.  Hence,  for  the  problem  at  hand  the 
necessary  conditions  are  also  sufficient. 

We  now  derive  the  differential  equations  for  the  adjoint  functions. 
Differentiating  with  respect  to  c  we  obtain 


:  ~«W‘. 
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4=  m 

*  dc 


-(♦  -  1)R 


or 

$  +  Rf  =  R 

Rt 

Multiplying  through  this  equation  by  the  integrating  factor,  e  ,  integrating 
from  t  to  T,  and  using  the  fact  that  4(T)  =  0,  we  obtain 

(25)  t(t)  =  1  -  eR(T”t) 

As  an  interpretation  of  th’s  function  we  note  that  the  expression 
i  -  ♦  «  ,R<T-*> 

is  the  future  (time  T)  value  of  a  dollar  gained  now  (timp  t). 

To  obtain  the  differential  equation  for  0  we  differentiate  with 
respect  to  w; 


(26)  0  «  -  =  (*  “  Dh’(w)  =  -h’(w)eR(T"t) 

Integrating  from  t  to  T  and  using  0(T)  =  0  yields 

(27)  0(t)  -  J**  h'(w)eR(T't)dt 

The  interpretation  of  0(t)  is  that  it  is  the  total  cost  per  unit  of  wheat 
of  storing  from  t  to  T  a  small  increment  of  wheat. 

Let  us  now  substitute  these  results  back  into  the  Hamiltonian  function 
and  give  a  new  interpretation  to  We  have 

(28)  -  =  [Rc  -  h(w)]eR(T‘t)  +  [P  -  peR(T_t)  -  0]m 

We  see  that  has  two  terms,  one  not  involving  the  control  variable.  The 


first  term  does  not  involve  the  control  variable  and  represents  the  future 
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(time  T)  value  of  the  interest.  Re  received  on  the  cash  account  less  the 
storage  charges  -h(w)  on  the  wheat  account  at  time  t.  The  second  term 
involves  the  switching  function 

(29)  s  =  P  -  -  0 

which  can  now  be  interpreted.  It  is  the  future  (time  T)  price  of  wheat,  P, 
less  the  future  value  pe  v  of  paying  now  (time  t)  for  wheat  at  price  p, 

and  less  the  total  storage  cost  0  from  time  t  to  T  of  the  additional 
wheat  bought  now.  The  optimal  policy  given  in  (24)  can  now  be  stated  as: 
buy  (sell)  the  maximum  amount  at  time  t  if  it  will  yield  a  net  profit  (loss) 
at  time  T  where  the  net  profit  or  loss  is  to  be  evaluated  taking  into 
account  all  possible  repercussions  of  the  purchase  or  the  sale  decision.  If 
this  optimal  policy  is  followed  then  the  second  term  in  (28)  gives  the  total 
addition  to  future  profits  resulting  from  the  optimal  purchase  (or  sale)  of 
wheat  now. 

We  can  now  give  a  brief  interpretation  of  that  includes  all  of  the 

foregoing.  The  function  is  the  future  (time  T)  value  of  the  resent 

(time  t)  cash  flow  and  consists  of  two  parts:  the  first  part  is  the  cash  flow 
due  to  previous  (before  time  t)  decisions,  and  the  second  part  is  the  cash 
flow  due  to  present  (time  t)  decisions. 
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4.  INVENTORY  HOLDING  COSTS 

It  may  be  interesting  to  see  how  the  function  0 ( t: ) ,  the  future  value 
of  inventory  holding  cost,  behaves  under  some  specified  forms  of  the  function 
h(w).  First,  let  us  consider  the  case  where 

(30)  h(w)  =  Sw^ 

Since  the  optimal  policy  given  in  (24)  is  bang-bang,  namely  either  to  bu,r  at 
the  maximum  rate  or  to  sell  at  the  maximum  rate, 

(31)  w  =  m  =  M 

for  any  t  between  two  adjacent  switching  points,  where  M  is  a  constant  M 
or  M  depending  upon  the  sign  of  the  switching  function.  Remembering  this, 
we  substitute  (30)  into  (27)  and  integrate  (27)  by  parts,  obtaining, 


(32) 


0 ( t )  =  2S  JT  weR(T"T)dT 

t 


=  2S 


)e 


R(T-t) 


T 


t 


2S  , 
—  (we 


R(T-t) 


-  w(T))  + 


2SM 

2 


(e 


R(T-t) 


1) 


for  any  t  in  the  interval  t*  t  £  T,  where  t*  is  the  last  switching 
point.  By  solving  for  w(t)  in  a  stepwise  manner,  we  can  then  search  for 
all  switching  points. 

To  show  another  example,  let  h(w)  be  defined  as  follows: 

(33)  h(w)  =  S  Jw | 

Then 


I 


4 


(34) 


h'  (w) 


S  egn  w 
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where  sgn  v  Is  1  if  w  >  0,  is  0  if  w  =  0,  and  is  -1  if  w  <  0.  Then, 
for  any  t  in  t**  'Sts:  T  where  t**  is  the  last  time  when  the  sign  of 
w  is  changed,  the  function  0(t)  is  derived  as 

(35)  0(c)  =  h' (w)  JT  eP(T*T)dT=  h'(w)(-  eR(T“T) jt 

t 

=  h'(w)  R(T-t)  _  ^  _  Ssgn  w  R(T-t)  _  ^ 

R  R 

Given  the  function  0(t),  the  switching  function  becomes: 

(36)  .  -  ?  -  (eR<T-t)  -  1) 

K 

This  may  be  written  as 


(37)  .  =  eR(T-t)  ((CMSJt  +  P)e-R<I-t>  -  .  pI 

-  eR<T-'>  (u-p) 

where  the  function  u  is  defined  as 

os)  u  =  (  2-asut  +  P]e-R(T-t)  -  5-ap-s 


Since  the  sign  of  s,  which  is  of  our  interest,  is  equal  to  the  sign  of  u-p 
for  R  #  0  and  finite  T  and  t,  we  may  draw  a  chart  for  u  for  given  values  of 
R,  S,  T,  P,  and  sgn  w,  and  compare  it  with  p  to  obtain  sgn  (u-p). 

Let  us  consider  the  case  where 

R  =  .01  (17.  per  month  compounded  continuously) 

S  =  .20  (dollars  per  unit  per  month) 

T  =  12  (month) 


p(T)=  10  (dollars) 
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1; 
1 1 


Q  ■> 

Then,  if  we  let  u  ,  u  ,  and  u  be  the  function  u  in  which  sgn  w  =  +1,  0, 
and  -1,  respectively,  they  may  be  written  as  follows: 


+ 

u 


o 


u 


u 


30e-w-t> 

10e-R(T-t) 


-  20 


+  20 


In  Figure  1,  three  curves  are  prepared  for  these  functions It  can  be 
easily  seen  from  (38)  that  for  a  larger  inventory  holding  cost  S,  the  spreads 
between  u+  and  u°  as  well  as  u°  and  u  become  wider  for  any  t  except  t  =  T, 
and  that  for  a  larger  interest  rate  R,  the  same  spreads  become  smaller  and 
all  three  curves  become  more  convex  starting  from  lower  points  at  t=0. 

The  curves  in  Figure  1  show  the  switching  criterion  in  the  following 
manner.  Let  us  suppose  that  the  sign  of  inventory  balance  does  not  change 
throughout  the  entire  period  of  o  £  t  £  T.  Then  u+  or  u  in  Figure  1  is 
used  to  compare  with  p.  If  p  is  belcw  u  (u+  or  u  depending  upon  sgn  w)  for 
any  time  interval  t^  s  t  s  tj,  this  means  that  sgn  (u-p)  (hence  sgn  s)  is 
positive.  Thus,  by  (24),  m  is  set  equ-.l  to  M  during  this  period,  where  M 
is  the  maximum  purchase  rate.  On  the  other  hand,  if  n  lies  abov°  u  for 
any  time  interval,  m  is  ~et  equal  to  M  during  this  period,  where  M  is  the 
maximum  sales  rate.  If  p  coincides  with  u  for  any  time  interval,  the  total 
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For  0  s  x  s  .12, 


e 


•x 


can  be  approximated  by 


e'X  *  1  -  . 95x 

with  errors  less  than  .0013  in  the  absolute  values.  The  curves  in  Figure  1 
are  drawn  based  on  this  approximation. 


16 


value  of  the  firm's  assets  at  time  T  Is  unaffected  whether  the  firm  buys  or 
sells  wheat  or  do  nothing  during  this  period,  although  the  amount  of  the 
inventory  balance  at  time  T  is  certainly  affected. 

Figure  2  shows  some  examples  of  p  together  with  indicated  decisions 

on  m. 

Now  let  us  consider  situations  where  the  sign  of  w  changes  during  the 
period  of  0  s  t  i  I.  Consider  the  case  (c)  in  Figure  2  where  p  =  7.5  for 
0  £  t  £  6,  Suppose  that  the  maximum  purchase  and  sales  rates  (M  and  M)  are 
set  equal  to  1  and  -1,  respectively,  and  that  the  firm  has  only  2  units  of 

wheat  at  t  =  0.  Then  at  t  =  2,  the  inventory  balance  becomes  0,  hence  we 

o  ■f*  o 

must  compare  p  with  u  rather  than  u  .  Since  p  lies  below  u  ,  we  must  buy 

(rather  than  sell)  wheat.  However,  as  soon  as  we  start  buying  wheat,  w 

becomes  positive;  hence  p  must  be  compared  with  u+  which  lies  below  p, 

indicating  "sell".  On  the  other  hand,  we  cannot  keep  selling,  since  as  soon 

as  w  becomes  negative,  p  must  be  compared  with  u  which  indicates  "buy". 

Thus,  once  it  is  reached,  v  =  0  is  an  equilibrium  point  for  any  price  between 

u+  £  p  £  u  at  any  time  t.  If  the  initial  inventory  is  0  and  p  lies  entirely 

between  u+  and  u  ,  there  will  be  no  purchases  or  sales  of  wheat  during  the 

entire  period.  This  situation  is  shown  in  Figure  3  by  an  example. 

The  analysis  as  elaborated  in  the  above,  may  be  difficult  to  obtain  if 
the  inventory  holding  cost  is  expressed  as  a  more  complicated  function  of  w, 
since  the  switching  function  s  may  become  a  function  of  w,  rather  than  Just 
the  sign  of  w.  However,  such  a  problem  can  be  solved  at  least  numerically  by 
means  of  using  an  analogue  computer  or  by  means  of  transforming  the  problem  into 
a  discrete  problem  as  elaborated  in  our  paper  [3]  and  the  solution  obtained  by 
means  of  a  digital  computer. 


5  6 

Figure  2 
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5.  TRANSACTION  AND  SPOILAGE  COSTS 

In  this  section  we  work  out  a  couple  of  minor  variations  in  the  model. 
First  we  assume  that  thei?  is  a  constant  proportion  of  the  wheat  that  spoils. 
The  rate  of  spoilage  is  given  by  a  positive  constant  A.  Second  we  assume 
that  a  tiansaction  of  size  m  incurs  a  unit  tv  reaction  cost  f(m),  where 
f  is  a  monotone  increasing  function  for  al1  values  of  m  with  f(0)  =  0. 

Later  on  we  shall  assume  that  f  is  unbounded  above  and  belcv,  and  3how  that 
we  may  drop  the  bounds  (17)  wher  this  is  true.  Further,  we  shall  assume 
throughout  this  section,  w  ^  0  for  all  t  in  the  range  considered. 

We  start  by  retaining  equations  (16)  and  (17).  The  rate  of  change  of 
cash  is 

(39)  c  ~  Rc  -  h(w)  -  [p  +  f(m)]m 

the  change  from  (18)  being  only  in  the  addition  of  the  transaction  cost  f(m)m. 
Note  that  the  transaction  cost  f(m)m  ^  0.  Similarly,  the  rate  of  change  of 
wheat  is  given  by 

(90)  w  =  a  -  Aw. 

the  only  change  from  (19)  is  in  the  spoilage  term,  -Aw. 

Substituting  these  into  x,  we  have 

(91)  x  =  £  +  Pw  =  Rc  -  h(w)  -  PAw  -  (p  -  f(m)  -  P]m 
From  these  and  the  same  adjoint  functions  as  before  we  have 

(42)  =  (1  -  i)  [He  -  h(w)J  +  (0  -  P)  Aw  +  [P  -  0  -  (l-t)(p+f  («))  !m 

The  differential  equation  for  t  is  exactly  the  same  as  before,  and  hence  j 
is  given  by  (25).  From  this  we  obtain  an  expression  for  that  is, 

(43)  -  =  (Rc  -  h(w)]eR^T  "  +  (0  V')kv  +  [P-0-eR^  C ^ (p+f (m) )  )m 


which  hau  similar  interpretations  as  before.  Namely,  the  first  term  in  (43) 
is  the  currant  change  in  cash  due  to  past  decisions;  the  second  term  is  the 
currant  lose  disc  to  spoilage;  the  last  term  iiwolves  the  switching  function 

(44)  a  ~  P  -  0  -  eR^  C/I(p  +  f(m)) 

which  is  just  like  (29)  with  transaction  costs  added  in.  Thus  the  decision 
ns  to  whether  or  not  to  buy  now  takes  into  account  the  transaction  cost.  The 
optimal  policy  is  again  bang-bang  and  is  described  in  (24),  provided  K  and  M 
are  not  too  large  in  absolute  value.  We  return  to  this  point  later. 

The  derivation  of  the  adjoint  function  0  is  a  little  more  complicated. 
It's  differential  equation  is 

(45)  0  -  A0  =  -  h'  (w)e^T-t^  -  PA 

Multiplying  by  the  integrating  factor  e  At,  integrating  from  t  to  T,  and 
using  the  boundary  condition  0 (T)  =  0  yields 

(46)  0 ( t )  =  eAt+RT  JT  h'(w)e"(A+R)TdT  +  P(1  -  e'A(T"t)) 

t 

The  first  term  Is  just  like  (27),  except  that  the  interest  rate  R  has  been 
augmented  by  the  spoilage  rate  A.  (In  this  model,  we  are,  in  effect,  being 
charged  for  storing  the  spoiled  wheat;  that  is  an  indirect  storage  cost  that 
appears  here  as  an  increase  in  interest  rate.)  The  second  term  of  (46) 
measures  the  actual  net  weight  of  wheat  taking  into  account  the  spoilage 
that  will  occur  between  time  t  and  T  if  we  have  a  bushel  of  wheat  now. 

Let  us  now  return  to  the  switching  function  (44)  and  the  fcrm  of  the 
optimal  policy.  Assume  that  f(m)  is  unbounded  above  and  below  and  is  mono¬ 
tone  increasing  (it  need  not  be  continuous),  and  also  f(0)  =  0.  Suppose 
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that  we  do  not  have  bounds  (17)  on  m.  In  Figure  4  we  plot  the  situation 

RfT-t) 

for  two  different  prices  p*  with  p*e  <  P  -  0  and  p**  with 

R(T-t) 

p*'*e'  >  P  -  0.  In  the  first  case,  it  clearly  pays  to  buy  wheat  and 


the  optimal  strategy  is  buy  m*  bushels  of  wheat  where  m*  solves  the 
equation 

(47)  (p  +  f(m))eR(T_t)  =  P  -  0 

In  the  second  case,  the  optimal  strategy  is  to  sell  m**  units  where  m**  is 
the  solution  of  (in  this  ca3e  the  negative  of)  (47).  The  fact  that  there  are 
solutions  of  (47)  in  each  case  follow  from  the  assumption  that  f(m)  is  un¬ 
bounded  above  and  below. 

In  the  case  in  which  (17)  holds  as  well  it  is  necessary  to  compare  the 
value  m*  with  M  and  choose  whichever  is  smaller  to  determine  the  optimal 


I 


i 
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purchase;  and  also  compare  m**  with  M  and  choose  whichever  is  larger 

(least  negative)  to  determine  the  optimal  sale.  j 

i 

i 

* 

6.  CONCLUSIONS  j 

We  have  oemonatrated  in  the  above  hcFw  mathematical  control  theory  tnav  be  J 

profitably  applied  to  the  overall  operations  of  a  firm  via  the  accounting  and  ! 

budgeting  system,  although  the  example  we  used  is  an  extremely  simple  one  in  ! 

order  to  obtain  the  theoretical  insight  in  the  clearest  manner.  It  is  cer¬ 
tainly  true  that  as  we  increase  the  number  of  variables  (or  accounts  in  the 
accounting  sense)  the  problem  becomes  harder  and  harder  to  solve.  It  is  also 
clear  that  for  practi  1  applications,  we  must  formulate  the  problem  in  a 

discrete  manner  in  order  that  it  may  be  solved  on  a  computer.  However,  we 

I 

are  leaving  these  problems  to  our  follow-up  paper  [3]  as  we  mentioned  at  the 

beginning.  In  the  present  paper,  we  are  primarily  interested  in  the  j 

I 

theoretical  clarification  of  the  problem.  j 

As  a  result  of  applying  the  mathematical  control  theory,  riot  only  can  J 

we  develop  a  schedule  for  optimal  control  in  terms  of  the  values  of  control  ‘ 

variables  but  also  we  can  prepare  optimal  projected  financial  statements 
(or  budgets)  similar  to  the  one  discussed  in  the  paper  by  Ijiri,  Levy,  and  * 

Lyon  [2 j .  Then  any  departure  of  the  actual  from  the  planned  activities  may 
be  analyzed  in  terms  of  errors  in  the  model  and  suboptimal  performances. 

There  is  one  additional  factor  derived  from  the  model  that  is  of 
significant  interest  from  the  accounting  standpoint.  That  is  the  accounting 
meaning  of  the  adjoint  variables.  There  have  been  several  attempts  in 
accounting  valuation  theory  to  use  the  discounted  future  cash-flow  as  a  basis 


of  asset  valuation.  As  we  have  seen  above,  the  adjoint  variables  offer  much 
richer  basis  of  valuation,  since  it  takes  into  account  not  only  interest  but 
Iso  the  effect  of  an  incremental  unit  of  an  asset  upon  the  objective  of  the 
firm  taking  into  account  the  entire  interaction  of  the  firm’s  activities.  We 
are,  of  course,  not  advocating  the  use  of  the  adjoint  variables  as  a  sub¬ 
stitute  for  all  valuation  basis;  however,  they  do  provide  useful  informa¬ 
tion  for  the  planning  purposes  of  the  management,  in  an  analogous  way  as  the 
dual  evaluation  in  linear  programming  does. 

For  instance,  in  the  case  discussed  in  Section  2,  the  adjoint  variable 
ijt  that  is  attached  to  the  variable  c  for  cash  provides  a  criterion  whether 
a  further  investment  in  this  firm  is  more  profitable  than,  say,  an  investment 
in  a  bond  which  will  yield  at  a  certain  rate.  A  similar  comparison  may  be 
made  by  the  use  of  the  adjoint  variable  for  wheat.  Furthermore,  the  switching 
function  s  may  be  used  as  a  criterion  for  expanding  the  upper  and  lower  bound 
for  wheat  transactions  since  the  terminal  asset  value  is  increased  by  |sj 
if  the  bound  is  extended  by  one  unit. 

We  are  far  from  having  exploited  this  fruitful  area  o^  applications. 

This  paper  is  only  the  beginning  of  many  interesting  and  useful  concepts  and 
techniques  that  may  be  derived  from  the  applications  of  the  mathematical  control 
t!  ry  in  accounting  and  budgeting. 


24 


BIBLIOGRAPHY 


[1]  Gelfond,  I.  M.  and  S.  V.  Fomin,  Calculus  of  Variations,  Englewood 
Cliffs,  New  Jersey,  Prentice-Hall,  1963. 


i2]  Ijiri,  Y.,  F.  K.  Levy,  R.  C.  Lyon,  "A  Linear  Programming  Model  for 

Budgeting  and  Financial  Planning,"  Journal  of  Accounting  Research, 
Vol.  1,  No.,  2,  Autumn  1963,  pp.  198-212. 


[3]  Ijiri,  Y.  and  G.  L.  Thompson,  "Applications  of  Mathematical  Control 
Theory  to  Accounting  and  Budgeting:  III.  The.  Discrete  Wheat 
Trading  and  Discrete  Joint  Production  Models,"  forthcoming. 


[4]  Pontryagln,  L.  S.,  R.  V.  Boltyanski,  R.  V.  Gamkraldze,  and  E.  F. 

Mishchenko,  The  Mathematical  Theory  of  Optimal  Processes, 
John  Wiley  and  Sons,  New  York,  1962. 


[5]  Tou,  J.  T.,  Modern  Control  Theory.  McGraw-Hill  Book  Co.,  New  York, 
1964. 


DOCUMENT  CONTROL  DATA  -  R  &  D 

(Security  classification  ot  title,  body  of  abstract  and  indexing  annotation  must  be  entered  when  the  overall  report  Is  classified) 


t  originating  activity  (Corporate  author) 


Graduate  School  of  Industrial  Administration  i2b.  croup 
Carnegie-Mellon  University  I 


3  REPORT  TITLE 


APPLICATIONS  OF  MATHEMATICAL  CONTROL  THEORY  TO  ACCOUNTING  AND  BUDGETING 
I,  THE  CONTINUOUS  WHEAT  TRADING  MODEL 


A  DESCRIPTIVE  NOTES  ('type  of  report  andjnclueive  dates) 


5  au  THORI3)  (first  name,  middle  initial,  last  name) 


Yuji  Ijiri 
Gerald  L.  Thompson 


8 a.  CONTRACT  OR  GRANT  NO. 


9a,  ORIGINATOR'S  REPORT  NUMBER!*) 


NONR  760(24) 


6.  PROJEC  T  NO 


Management  Sciences  Research  Report  No.  152 


NR  047-048 


9b.  OTHER  REPORT  NO(S)  (Any  other  number  a  that  may  be  me  signed 
thie  ii sport) 


10  DISTRIBUTION  STATEMENT 


Distribution  of  this  document  is  unlimited 


I  2  SPONSORING  MIL*  TARV  ACTIVITY 


Logistics  and  Mathematical  Statistics  Br. 
Office  of  Naval  Research 
Washington,  D.  C.  20360 


13.  ABSTRAC  T 


A  brief  introduction  to  continuous  mathematical  control  theory  is 
presented  £4*#t.  Then  a  model  having  two  state  variable  accounts,  cash  and 
wheat,  is  defined  by  me  ns  of  differential  equations.  Adjoint  functions, 
the  Hamiltonian  function,  and  the  optimum  (bang-bang)  policy  are  derived. 
Interpretations  of  the  Hamiltonian  and  switching  functions  and  a  numerical 
example  are  then  given.  Finally  the  results  are  extended  to  the  case  where 
transaction  costs  and  spoilage  costs  are  incurred. 


DD  ,?o?..1473  < 


PAGE  I  ) 


Unclassil 


S/N  0101 -807.68) I 


ecurity  ficAhon 


Unclassified 

Security  fiiitificiUon  * 

Control  Theory 
Maximum  Principle 
Budgeting 
Accounting 

Wheat  Trading  Example 


DD  .'•T..1473  <back> 


